# Full text of "LHC physics: challenges for QCD"

## See other formats

LHC PHYSICS: CHALLENGES ^ FORQCD o o . ^ ■ Stefano Frixione > o ^ . INFN, Sezione di Geneva, Via Dodecanese 33, 16146 Geneva, Italy OO ' E-mail: Stefano. Frixione© cern.ch • I review the status of the comparisons between a few measurements at hadronic , colliders and perturbative QCD predictions, which emphasize the need for improv- ed ' ing the current computations. Such improvements will be mandatory for a satis- factory understanding of high-energy collisions at the LHC ■ One of the main goals of the LHC will be the search of the Higgs bo- ^ . son, the only piece of the otherwise thoroughly tested Standard Model (SM) QhI of which we miss experimental evidence. Regardless of the existence of Q^' the Higgs boson, LHC will likely shed light on the physics Beyond-the-SM 1^ . (BSM), which should be within reach, as LEP results and solar and atmo- '"^ ' spheric neutrino data suggest. In this context, the role played by QCD may appear an ancillary one. However, this seems to be quite at odd with the fact ^ ■ that LHC is a hadronic collider, and strong interactions will be responsible for a prominent part of the reactions taking place. On the other hand, one may claim that an accurate knowledge of QCD predictions will not be nec- essary for the discovery of - say ~ the Higgs in presence of a striking feature such as a narrow mass peak, which could give the possibility of normalizing the background directly with the data. However, if the discovery is based on a counting experiment, it relies by definition on precise predictions for SM backgrounds, dominated by QCD. These predictions will also prove essential for after-discovery studies, when the properties of the newly-found particles will have to be determined, or in the case of absence of BSM signals, in order to set limits on BSM scenarios. Having argued that QCD studies are one of the keys for a successful LHC physics program, one question remains: are there still motivations to keep 1 on working on them? The answer is again yes. The question is legitimate, since the striking success of perturbative QCD in predicting experimental results may lead to think that the current knowledge is sufficient to tackle the problems posed by the LHC. Unfortunately, this is not the case: be- cause of the large energy available, hard reactions will occur in a previously unexplored regime, with peculiar kinematic characteristics, such as the si- multaneous presence of many, well-separated hard jets. These new features require to improve QCD predictions, either by increasing the accuracy of the computations, or by considering issues which could have been safely ne- glected up to now. Although one could give theoretical arguments for the improvements that need to be achieved, it is also instructive to look at those measurements which so far could not be described by perturbative QCD in an entirely satisfactory way, which in fact, rather than hinting to a funda- mental problem of QCD, indicate that a deeper understanding is desirable of some aspects of the computation. This is in fact relevant to the problem of QCD predictions for the LHC, since aspects which may be marginal in the current phenomenological picture will be much more important in the future. loO ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 5 10 15 20 25 Pt (GeV) Figure 1: data versus theoretical predictions Let me start with what has been seen historically as a major problem 2 in QCD, namely single-inclusive b production at hadron colliders: although NLO QCD appears to give a good description of data in terms of shape, the rate is typically underestimated by a factor of 2. Recently, this picture seemed to receive strong support from a CDF measurement PP, where the average of the data/theory ratio for mesons was quoted to be 2.9 it 0.2 it 0.4. However, it has been subsequently pointed out [2] that this value is mainly due to an improper treatment of the theoretical prediction, and that the correct result is 1.7 it 0.5 it 0.5. This value originates from using the state-of-the-art computation for single- inclusive b Pt spectrum (FONLL and especially from the observation that the fragmentation function b ^ B as presently determined by using e^e~ data is not particularly appropriate for the case of hadronic collisions. It is interesting to note that the proper choice of the fragmentation parameters leads to a comparison between theory and data (see fig. ^ which is basically identical to that relevant to 6-jet transverse energy [^[3], i.e. to an observable independent of the details of the fragmentation mechanism. Therefore, b production at hadron colliders should not be regarded any longer as a reason of concern; on the other hand, the size of the theoretical uncertainties (solid band in fig. ^ prevents more stringent tests. It is likely that a major source of improvement in this respect would be the computation of the b cross section to NNLO accuracy; a further enhancement of the rate should be expected from small-x jOl EI and threshold resummations. Finally, let me remind that results for total rates for b production in ep, and 77 collisions are in such a disagreement with NLO QCD predictions, that no viable explanation for this discrepancy has been found in any BSM scenario. It is necessary to note that in many cases the experimental results are extrapolated to the full phase space from a rather narrow visible region. However, it is encouraging that in a few cases the data are also presented without the extrapolation outside the visible region, and in this way they are fully compatible with QCD predictions. This points out that some problem may be hidden in the Monte Carlo (MC) simulation of heavy flavour production, especially in the low region. Let me now turn to jet production at the Tevatron. The excess of CDF single-inclusive jet data [Sj over NLO predictions at large has been re- garded with much interest, being a potential signal of new physics. On the other hand, it has been immediately observed that, by suitably adjust- ing the gluon density in the proton, one can obtain sets of PDFs (denoted with the "HJ" suffix by CTEQ) which give a decent global fit, and result in predictions for jets compatible with CDF measurements. This proves that QCD has enough flexibility to accommodate the excess, but it is disturbing that the "HJ" family is not the preferred one according to a (unweighted) 3 global fit procedure. This situation changes when DO jet data [TU] are in- cluded in the global fit, since the resulting PDF set turns out to belong to the "HJ" family (although it is not called accordingly). This happens be- cause the Bjorken x's relevant to high-px jet production are also relevant to 0.8 0.4 0.0 -0.4 5HJ 6M 5M 0<T|<0.5 1.0 0.5 0.0 -0.5 DO data 0<T1<0.5 — PDF uncertainty 100 200 300 400 500 100 200 300 400 500 Figure 2: Data/theory ratios for single-inclusive jet production. From ref. small-px; large- r/ production, a region probed by DO (see ref. for a dis- cussion). The comparison between theory and data is now fully satisfactory: the left panel of fig. El presents the ratio of the DO data ^10, over theoretical predictions obtained with the PDF set CTEQ6M1 for reference, the ratios of theoretical predictions obtained with other PDFs are also given; the situation for CDF data is analogous. Although this excludes any evi- dence of new physics in this channel, the situation may still change. This is shown in the right panel of fig. |21 where the band represents the span of the theoretical predictions due to PDF uncertainties. Thus, single-inclusive jet measurements at the Tevatron clearly stress the importance of precise PDF determinations (which implies the necessity of using data from many comple- mentary processes), and of accurate estimates of the uncertainties affecting the PDFs. It is also worth reminding that there are at least a couple of unpleasant aspects of jet production at colliders. Firstly, by reconstructing jets using a fcx-algorithm, DO ^2] finds large discrepancies with QCD predic- tions, a fact difficult to reconcile with the observation that, at the NLO and for suitable choices of jet-recombination parameters, no major differences can be seen in the single-inclusive px distribution of jets reconstructed with the /cx or the cone algorithms. Although more experimental analyses with the kx algorithm are necessary in order to confirm the result of ref. jl2j . it is interesting to observe that hadronization corrections and the underly- ing event modelling, as simulated by MC's, affect fairly differently the jets reconstructed with different algorithms. Secondly, the ratio of jet cross sec- tions measured by CDF and DO at different cm. energies {V^ = 1.8 and 4 0.63 TeV) and fixed = 2pT,/^/S, is not in agreement with QCD (which can predict this quantity in a fairly accurate manner); besides, CDF and DO data are also mutually incompatible at low Xt- However, the differences are largely within the uncertainties due to neglected power-suppressed effects. In both cases, it appears that the experimental analyses would benefit from a much deeper understanding of the interplay between perturbative and non-perturbative physics, both at the level of MC's (whose use to estimate the hadronization corrections applied to NLO predictions is rather empiric), and in the context of approaches analogous to that of DMW |13l I14| . O 1.5 0) o _^0.5 \ 2 o Q ^0.5 \ -1 ▲ CDF 1800 GeV Data • CDF 630 GeV Data 10 10 Photon Xt |7ll<0.9 NLO QCD (CTEQ5M) \s.=E^ 1.6<|r||<2.5 ■ DO daUi 0.02 0.04 0.06 0.0 0.1 0.12 0.14 Figure 3: Comparison between prompt-photon data and QCD predictions Let me finally deal with the case of 7 production. One of the reasons of interest in this process is that, through the leading order partonic reaction qg 'jq, it can in principle constrain the gluon density in a rather clean way. This is even more interesting if one observes that the typical range in probed by fixed-target experiments is quite similar to the range in x^^^^ probed in moderate- and large-pT jet production at the Tevatron. Unfortu- nately, as discussed in ref. fixed-target photon data are not mutually compatible. This nowadays prevents the use of these data in global PDF fits. The problem is not peculiar to fixed-target experiments, since the agreement between isolated-photon data at the Tevatron and QCD predictions has al- ways been pretty marginal, as shown in the left panel of fig. Ol where CDF data are compared to NLO QCD results ^Zj. The situation, partly be- cause of larger statistical errors, improves for DO ^Sl) as shown in the right panel of fig. 01 (here, ratios of cross sections, rather than cross sections, are presented). It must be stressed that the isolation cuts used by experiments 5 at colliders are unlikely to be fully equivalent to those applied in theoretical computations. The whole situation is rather inconclusive; if QCD predic- tions have to be used in the context of analyses involving photons, it is urgent to devise a strategy that allows one to apply the same cuts at the theoretical and experimental levels. For example, in a high-energy environ- ment experiments prefer to define isolated photons with narrow cones; on the other hand, it is known that the isolated-photon cross section at NLO becomes unphysical (being larger than the fully inclusive one) for small iso- lation cone sizes U^l- The computation of yet higher orders, and especially of the all-order resummation of the relevant terms, would be necessary in order to solve this problem. Although no major problem for QCD emerges from the phenomeno- logical considerations given above, the picture would certainly look firmer (not necessarily better) if we had: a) NNLO results"*^; b) resummed results, matched with fixed-order ones^; c) better determinations of PDFs and their uncertainties; d) better understanding of the interplay between perturbative and non-perturbative physics; e) better models for the underlying event in MC's. It should be clear that these issues are of primary importance for LHC physics: large numbers of hard jets (and large K- factors), strong im- pact of the accompanying non-hard events, and many-scale processes will occur plentiful. This is also worrying in view of the fact that MC's, the ubiquitous tools in experimental analyses, are not reliable when an accurate description of large-angle emission is needed. It seems therefore mandatory to add to the list above: f) better MC simulation of hard emissions. Items a)-f) are all tough problems, but it is generally believed, thanks to recent developments, that substantial progress will be made before LHC comes into operation. In view of its relevance to experimental collaborations, in the following I'll concentrate on item f). It is useful to briefly remind how an MC works: for a given process, which at the LO receives contribution from 2 ^ no reactions, (2 + no)-particle configurations are generated, ac- cording to exact tree-level matrix element (ME) computations. The quarks and gluons (partons henceforth) among these primary particles are then al- lowed to emit more quarks and gluons, which are obtained from a parton shower or dipole cascade approximation to QCD dynamics. To lessen the impact of this approximation on physical observables, one can devise two strategies. The first aims at having extra hard partons in the final state; thus, in the example given above, the number of final-state hard particles would increase from no to no + ng. This approach is usually referred to as ^For certain processes and observables. 6 matrix element corrections, since the MC must use the (2 + no + raE)-particle ME's to generate the correct hard kinematics. The second strategy also aims at simulating the production of no + hard particles, but improves the computation of rates as well, to N^^LO accuracy. I'll generally denote the resulting MC as N"eLOwPS. There are basically two major problems in the implementation of ME corrections. The first problem is that of achieving a fast computation of the ME's themselves for the largest possible no+ng, and an efficient phase-space generation. A variety of solutions exist nowadays for this problem, imple- mented in packages which I'll denote as ME generators. The second problem stems from the fact that multi-parton ME's are IR divergent. Clearly, in hard-particle configurations IR divergences don't appear; however, the def- inition of what hard means is, to a large extent, arbitrary. In practice, hardness is achieved by imposing some cuts on suitable partonic variables, such as Pt's and (77, (/3)-distances. I collectively denote these cuts by 6sep- One assumes that n hard partons will result (after the shower) into n jets; but, with a probability depending on 6sep, a given n-jet event could also result from n + m hard partons. This means that, when generating events at a fixed no + number of primary particles, physical observables in gen- eral depend upon 5s£p] I refer to this as the J^ep'bias problem. Any solution to the Jsep-bias problem implies a procedure to combine consistently ME's with different no + n^'s. It should be stressed that, in presence of a 5sep bias, the interface of an ME generator (which is responsible for producing the hard configurations, i.e. the initial conditions for the shower), and a parton shower code is not, strictly speaking, an event generator, since the events depend somehow on the value of 6sep- In practice, the dependence is of the order of 20%, which is acceptable if one considers that, without ME corrections, multi-jet configurations predicted by standard MC's are com- pletely unreliable. A solution to the Jsep-bias problem has been presented, for e'^e~ collisions, in ref. ^20"^ (CKKW henceforth), and subsequently ex- tended (without formal proof) to hadronic collisions in ref. Loosely speaking, CKKW achieve the following: if an n-jet observable is affected by the 5sep bias in the following way k by applying the CKKW prescription one gets «r'(5.%+E^fc«siog2'"''^-p). (2) k 7 The implementation of CKKW in popular event generators for hadronic collisions in under way, and it will have reached a mature and well-tested stage when LHC will come into operation. The implementation of an N"^LOwPS can be seen as an upgrade of ME corrections: not only one wants to describe the kinematics of uq + hard particles correctly, but the information on N^^LO rates must also be in- cluded. First attempts at solving this problem have only recently become available, and only for the case Ue = 1. The striking feature of an NLOwPS is the computation of loop diagrams (which are necessary in order to com- pute total rates to NLO accuracy); this in general implies the presence of negative weights. This is a new feature in MC's, which however doesn't spoil their probabilistic nature. In fact, in NLOwPS the distributions of positive and negative weights are separately finite, at variance with what happens in NLO computations; thus, each of them can be unweighted and evolved separately, since no cancellation between large numbers is involved in this procedure. On the other hand, the contribution of loop diagrams im- plies that the 5sep-bias problem which affects ME corrections is simply not present. At the moment, the following codes implement different prescrip- tions for NLOwPS in hadronic collisions: $-veto MCONLO |2S1I211, GRACE_LLsub $-veto is based on the slicing method, and features Z* and W* production; it is affected by double counting according to the definition of ref. j23j . but numerically this problem seems to be of minor importance; it is interfaced with Herwig and Pythia. MC@NLO is based on the subtraction method, and features W~^W~, ZZ, WZ, bb, tt, SM Higgs, Z, W, and 7 production; it is interfaced with Herwig. GRACE_LLsub is based on the slicing method and on the fully-numerical computation of the matrix elements; it features Drell-Yan production, and is not affected by double counting only if the parton shower of ref. [25| is adopted (i.e., other showering codes cannot be used at the moment). The field of NLOwPS, still behind that of ME corrections, is rapidly evolving, and more ideas will appear in the future; soon, more processes will be implemented, and a thor- ough comparison between the various approaches will have to be made. In summary, a lot of interesting developments are currently occurring in QCD, which will provide a solid benchmark for LHC studies. It will be vital for experimental collaborations to exploit these results, both by using new MC tools (with ME corrections and NLOwPS) in the course of their analyses, and by considering the most precise theoretical results available, at fixed-order (NNLO) or in resummed computations (with NLL or NNLL accuracy). In the coming years, it will also be crucial to learn from the experience of HERA and the Tevatron, which will hopefully provide the 8 necessary data to determine PDFs at an unprecedented level of accuracy, and to test the various models for underlying events. It's a pleasure to thank the organizers for an interesting meeting; as a phenomenologist, I think it is beneficial for both theoretical and experimen- tal communities to have frequent (elastic) interactions. References [1] D. Acosta et al. [CDF C ollaboration], Phys. Rev. D 65 (2002) 052005 ja?Xw:hep-ph/01li359l . [2] M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 122003 |arXiv:hep-ph /0204025] . [3] M. Cacciari , M. Greco and P. Nason, JHEP 9805 (1998) 007 |arXiv:hep-ph /9803400] . [4] S. Frixione and M. L . Mangano, Nucl. Phys. B 483 (1997) 321 |alxi7:hep-ph/9605"270l . [5] B. Abbott et al. [DO Collaboration], Phys. Rev. Lett. 85 (2000) 5068 |arXiv:hep-ex /000802l] . [6] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366 (1991) 135. [7] J. C. Collins and R. K. Ellis, Nucl. Phys. B 360 (1991) 3. [8] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 77 (1996) 438 rarXiv:hep-ex/9601008 . [9] J. Huston, E. Kovacs, S. Kuhlmann, H. L. Lai, J. F. Owens, D. E. Soper and W. K. Tung, Phys. Rev. Lett. 77 (1996) 444 |arXiv^hep^ph/95 1 1386] . [10] B. Abbott et al. [DO Collaboration], Phys. Rev. D 64 (2001) 032003 |arXiv:hep-ex/0012046] . [11] D. Stump, J. Huston, J. Pumplin, W. K. Tung, H. L. Lai, S. Kuhlmann and J. F. Owens, JHEP 0310 (2003) 046 | arXiv:h ep-ph/0303013|. [12] V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 525 (2002) 211 |aIXiv:hep-ex/010904T] . 9 [13] Y. L. Do kshitzer and B. R. Webber, Phys. Lett. B 352 (1995) 451 |arXiv:hep -ph/95 04219| . [14] Y. L. Dokshitzer, G. Marchesini an d B. R. Webber, Nucl. Phys. B 469 (1996) 93 arXiv:hep-ph/ 9512336| . [15] P. Aurenche, M. Fontannaz, J. P. Guillet, B. A. Kniehl, E. Pilon and M. Werlen, Eur. Phys. J. C 9 (1999) 107 1 arXiv:hep-ph/9811382| . [16] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 65 (2002) 112003 |arXiv:hep-ex/ 0201004] . [17] M. Cluck, L. E. Gordon, E. Reya and W. Vogelsang, Phys. Rev. Lett. 73 (1994) 388. [18] V. M. Abazov et al. [DO Colla boration], Phys. Rev. Lett. 87 (2001) 251805 |arXiv:hep-ex /0106026|. [19] S. Catani, M. Fontannaz, J. P. Guillet and E. Pilon, JHEP 0205 (2002) 028 |arXiv:hep-ph/ 0204023| . [20] S. Catani , F. Krauss, R. K uhn and B. R. Webber, JHEP 0111 (2001) 063 arX iv:hep-ph /01092311. [21] F. Krauss, JHEP 0208 (2002) 015 |arXiv:hep-ph7 0205283l . [22] M. Dobbs, Phys. Rev. D 65 (2002) 094011 |arXiv:hep-ph701 11234] . [23] S. Frixione and B. R. Webber, JHEP 0206 (2002) 029 arXiv :hep-ph/0204244 . [24] S. Frixione, P. Nason and B. R. Webber, JHEP 0308 (2003) 007 arXiv:hep -ph/03 05252 . [25] Y. Kurihara, J. Fujimoto, T. Ishikawa, K. Kato, S. Kawabata, T. Munehisa and H. Tanaka, Nucl. Phys. B 654 (2003) 301 arXiv:hep-ph/02122T6l . 10